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ABSTRACT 

A  pilot  engaged  in  a  strike  mission  against  targets  of  opportunity 
is  continually  confronted  with  the  decision  as  to  how  many  and  which  of 
his  available  weapons  should  be  expended  against  a  variety  of  targets 
which  are  encountered .  Associated  with  this  decision  problem  is  the 
problem  of  determining  an  optimal  load  for  such  a  mission  subject  to 
the  constraints  of  available  payload  and  number  of  weapons  stations  on 
the  aircraft.  Certain  assumptions  are  made  concerning  the  distribution 
of  targets  within  the  target  area  which  lead  to  a  dynamic  programming 
formulation  of  the  decision  problem.  This  yields  a  system  of  ordinary 
differential  equations  which  are  solvable  recursively.   In  addition  to 
the  dynamic  programming  model ,  a  sub-optimum  determination  of  a  "best 
load"  is  discussed.  Although  not  as  complete  or  precise  as  the  dynamic 
programming  method,  this  formulation  is  more  readily  adaptable  to  squa- 
dron level  decision  making. 
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1.   INTRODUCTION 

The  problem  of  selecting  weapons  for  an  aircraft  strike  mission 
and  using  those  weapons  effectively  has  been  documented  quite  exten- 
sively.  However,  present  publications  are  devoted  almost  exclusively 
to  the  situation  where  the  target  is  assumed  known  prior  to  the 
mission  and  the  aircraft  can  be  loaded  for  taat  target,  or  where 
several  different  aircraft  can  be  loaded  differently  to  cover  an 
unknown  spectrum  of  targets.   These  approaches  are  inadequate  for 
determining  weapons  requirements  against  targets  of  opportunity.   A 
pilot  engaged  in  such  a  mission  is  confronted  with  the  decision  as 
to  which  and  how  many  of  his  available  weapons  to  expend  when  any 
target  is  encountered.   When  various  types  of  targets  might  be  en- 
countered, each  witn  a  different  strategic  or  tactical  value,  he 
wishes  to  maximize  the  value  of  the  targets  he  can  expect  to  destroy. 
Associated  with  this  is  the  problem  involving  the  choice  of  weapons 
that  should  be  initially  loaded  onto  the  aircraft  for  such  a  mission. 
Here  the  problem  is  constrained  by  maximum  payload,  weapons,  and  wea- 
pons stations  availability.   It  is  desired  to  obtain,  within  these 
constraints,  (1)  the  optimum  feasible  weapons  mix  to  be  loaded  for 
use  against  an  uncertain  spectrum  of  targets,  and  (2)  a  doctrine  for 
the  expenditure  of  these  weapons  when  targets  are  encountered. 

The  model  developed  herein  attacks  these  problems  in  reverse 
order.  First,  an  optimum  firing  doctrine  is  generated  for  all  feas- 
ible loads.  Then  that  load  which  yields  the  maximum  expected  return, 


NAVWEPS  Report  8664,  "Preferred  Aircraft  Loadings  Based  on  the 
Predicted  Combat  Effectiveness  of  Conventional  Air-to-Surface  Weapons," 
NOTS  TP  3710,  1965,  (SECRET). 
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when  expended  in  accordance  with  the  optimum  doctrine,  becomes  the 
optimal  load. 

ASSUMPTIONS  IN  THE  MODELS 

Certain  assumptions,  necessary  to  the  model  development  are  made 
concerning  target  distributions  within  the  target  area  and  tactics  of 
the  strike  aircraft. 

It  is  assumed  that  there  are  a  finite  number  of  different  target 
types,  each  randomly  and  independently  distributed  within  the  target 
area  with  known  density.  For  example,  from  photo -reconnaissance  or 
other  sources  it  might  be  known  that  the  density  of  trucks  within  the 
target  area  is  x  per  square  mile.  The  number  of  targets  of  each  type 
is  assumed  to  be  large  enough  so  that  the  probability  of  future  en- 
counters is  unchanged  by  consideration  of  the  number  already  met. 
Further,  these  targets  are  assumed  to  be  of  a  unitary  nature,  charac- 
terized by  the  existence  of  single  shot  kill  probabilities,  as  dis- 
tinguished from  area  targets  characterized  by  probabilities  of  des- 
troying a  certain  percentage  of  the  area  or  value  of  the  target. 

The  assumption  is  made  that  a  ratio  scale  of  target  values  can 
be  specified.  For  example,  with  the  least  valuable  target  as  a  base 
with  value  V,  the  achievement  of  a  kill  on  a  target  type  i  is  worth 
7V,  target  type  j  is  worth  8V,  etc. 

For  defensive  reasons  the  strike  aircraft  is  restricted  to  a 
single  attack  on  any  target  which  might  be  encountered.  The  pilot 
may  not  drop  a  bomb  on  the  target,  assess  damage,  and  then  make 
another  attack  if  a  kill  has  not  been  achieved.  He  may  however, 
expend  any  number  of  weapons  on  this  single  pass  provided  they  are  all 
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the  same  type.   That  i6,  while  he  may  expend  three  500  lb.  bombs,  he  may 
not  drop  one  500  lb.  bomb  and  one  1,000  lb.  bomb.   This  restriction 
is  not  necessary  to  the  development  of  the  models,  but  rather  recog- 
nizes that  differences  in  aim  points  between  different  weapons  make 
this  mixed  drop  undesirable. 

METHOD  OF  SOLUTION 

Using  data  on  weapons  weights,  constraints  of  aircraft  payload 
capacity,  and  the  number  of  available  weapons  stations  on  the  plane, 
all  feasible  loadings  can  be  enumerated.   Initial  loads  are  then  con- 
sidered which  are  not  subsets  of  some  other  feasible  load.   That  is, 
if  five  weapons  of  type  one  and  three  weapons  of  type  two  constitute 
a  feasible  load,  then  four  weapons  of  type  one  and  three  of  type  two 
is  also  feasible  but  a  subset  of  the  former.   The  kill  potential  of 
some  subset  of  a  feasible  load  is  clearly  less  than  that  of  the 
original  load  and  is  therefore  not  considered  in  the  determination 
of  an  optimal  initial  loading. 

A  functional  equation  for  the  expected  value  of  targets  killed 
is  then  set  up  as  a  dynamic  programming  process  yielding  a  value  for 
kill  potential  for  each  of  the  feasible  loads  being  considered.   The 
maximum  kill  potential  from  amongst  these  loadings  determines  the 
optimal  load  for  the  strike  aircraft.   Additionally  in  the  solution 
of  this  dynamic  programming  process,  the  optimal  firing  strategy  for 
the  pilot  is  obtained  for  any  initial  load.  The  optimized  sequential 
decision  process  yields  a  policy  which  determines  as  a  function  of 
the  operating  time  and  weapons  remaining,  the  best  weapons  to  be  ex- 
pended in  an  encounter  with  any  of  the  various  target  types. 
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Dynamic  programming  provides  a  complete  and  very  precise  means 
of  determining  the  optimum  load  for  a  "targets  of  opportunity"  type 
mission.  Except  in  the  very  simplest  cases  however,  the  problem  is 
quite  complex  and  will  require  a  computer  for  solution.  In  the  majority 
of  operational  situations  such  equipment  will  not  be  available.  Thus 
an  approximate  method,  suitable  for  hand  calculation,  is  needed  for 
determining  the  initial  load. 

In  order  to  achieve  simplicity,  this  "hand  model"  must  sacrifice 
detail  and  thus  completeness.   In  this  respect  it  is  a  sub-optimal 
method.  The  simplified  method  is  best  suited  for  comparing  a  few, 
rather  than  all,  feasible  loads  for  use  against  a  selected  target 
list.  By  choosing  the  best  weapon  for  each  target  type  and  comparing 
the  expected  returns  for  various  mixes  of  these  weapons,  a  "best  load" 
can  be  generated. 
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2.      GENERAL  DISCUSSION  AND  DEVELOPMENT   OF  THE  MODELS 

The   random  distribution  of   targets    leads    to  an  exponential   inter- 
encounter   time  distribution   for   the   strike  aircraft   searching   the 
target  area.        (All   targets  which  are  within   the  visibility  range  of 
the  aircraft's   track  are  assumed   sighted.      Alternatively   the   sightings 
could   be   reduced   by  some   uc„.aaation   factor  which   takes    into  account 
the  difficulty  associated  with  perceiving  a   given   target   type.      This 
however  would   not   alter   the  analysis  which   follows.)      The  mean  rate 
of   sightings    (the   reciprocal   of    the  mean  of   the   exponential  distri- 
bution of    inter-encounter    times)    is   denoted   by   r. 

DYNAMIC    PROGRAMMING  MODEL 
The  development   of   the  dynamic  programming   formulation  of   the 
problem   is   best   seen  by   first   considering   the   simpler   cases.      To   begin 

with,    a  model  has   previously  been  developed    for   the  allocation  of  wea- 

2 
pons  over  a  limited  time  interval.  This  can  be  utilized  for  gener- 
ating the  optimal  firing  strategy  in  the  simplest  case  of  one  target 
type-one  weapon  type.  The  optimal  loading  problem  is  not  pertinent; 
for  a  single  weapon  type  the  best  load  is  clearly  the  maximum  number 
that   can  be   carried   on   the  aircraft.      First  define: 


f   (t)    =  expected   number  of   targets   killed   in   the   time 
n  interval   (0,t)    if   there  are  n  weapons  at   time 

t  =  0  and   these  are  utilized  optimally. 


1 
Appendix  I 

2 

Bram,  J.,  "Allocation  of  Weapons  to  Targets  with  Exponential 
Arrival  Times  in  a  Limited  Time  Interval,''  OEG  Interim  Research 
Memorandum  32,  1963. 
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An  initial   loading  of  N  weapons  and  a   time  T  allowable   for   search 

is  assumed,   and   f   (t)   will  be  calculated   for  0  ^  n  ^  N  and  0  ^  t  ^  T. 
n 

Define: 


p  ■  probability  that  j  weapons  expended  against  a  target 
j   will  achieve  a  kill. 


J 


,,1 


Clearly  f^(t)  is  dependent  on  the  p  '  s  and  on  r,  the  encounter 
rate. 

By  Bellman's   "principle  of  optima  lity,"J'  for  0  £  t-^  t,    the 
optimal  policy   in  the   interval   (0,    t)  ,   given  n  weapons  at  time   t  =  0, 
must  also  be  optimal   in  the  subinterval   (t, ,    t)    given  n  -   n     weapons 
at   time  t  =  t    ,  whare  n     weapons  were  expended   in  the   first  subin- 
terval  (0,    t.). 

Let  t  =  At  «  T.  By  definition  of  the  exponential  arrival  pro- 
cess, the  probability  of  an  encounter  in  a  small  interval  of  time  At 
is  proportional  to  At  with  constant  of  proportionality  r,  the  en- 
counter rate.  Thus  the  probability  of  one  encounter  in  the  time 
interval  (0,  At)  is  rAt  and  no  encounters  occur  with  probability 
1  -  rAt. 

Breaking  the   time   interval   (0,    t)    into  the  two  sub-intervals 
(0,   At)    and   (At,    t)    yields: 


(1  -  rAt)  =  the  probability  that  no  encounters  will  occur 
in  the  first  sub- interval. 

p  =  the  expected  return  in  the  first  sub-interval 
-*   having  expended  j  weapons. 


1 
Bellman,  R.,  Dynamic  Programming,  (Princeton,  New  Jersey: 

Princeton  University  Press,  1957),  p.  83. 
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f  (t  -  At)  »  the  expected  return  in  the  second  sub- interval 

having  expended  no  weapons  in  the  first  and  following 
an  optimal  strategy  in  (At,  t)  . 

f   (t  -  At)  -  expected  return  in  the  second  time  interval  having 
n-j  expended  j  of  the  n  weapons  in  the  first  sub- 

interval  and  behaving  optimally  in  (At,  t) . 


The  expected  return  over  the  whole  interval  (0,  t)  is  then: 

(1  -  rAt)  f  (t  -  At)  +  rAt  (p  If    (t  -  At)) 
n  j    tt-j 

To  optimize  over  the  period  (0,  t) ,  j  must  be  chosen  such  that  this 
expression  is  a  maximum.  Therefore 


(1)      f  (t)  =  (1  -  rAt)  f  (t  -  At)  +  rAt  max   (p  +  f    (t  -  At)) 
n  n  i    n-j 

Uj*n 


The  maximization  here  is  over  the  range  l^j^n  rather  than  O^j^n 
since  it  clearly  would  never  pay  to  let  a  target  go  by  without  firing 
any  weapons,  saving  them  all  for  future  encounters  which  are  not 
certain  to  occur. 

Rearranging  (1),  dividing  by  At  and  taking  the  limit  as  At 
approaches  zero  yields: 


(2)       f'(t)  -  -  rf  (t)  +  r  max   (p  +  f    (t)) 
n  n  n- j 

From  (2),  f  (t)  can  be  calculated  for  all  n  and  t  recursively 
n 

as  follows:   By  definition,  f  (t)  =  0  for  all  t.   (The  expected 

0 
return  is  zero  over  any  period  of  time  if  there  are  no  weapons  with 

which  to  operate.)  Also  fn(0)=  0  for  all  n.   (The  expected  return 

is  zero  if  there  is  no  time  in  which  to  operate.)   Solving  (2)  for 
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n  =  1  yields: 

fx(t)  -  px  (1  -  e"") 

Note  that  this  result  was  to  be  expected.   Exponential  inter- encounter 
times  lead  to  a  Poisson  counting  process  for  the  number  of  encounters 
within  a  time  interval.  The  term  (1  -  e   ),  is  simply  the  prob- 
ability of  at  least  one  encounter  in  the  time  interval  (0,  t)  ,  and  p 
is  the  probability  that,  given  an  encounter,  a  kill  will  be  achieved 
with  the  expenditure  of  the  single  weapon. 

Utilizing  f  (t)  equation  (2)  can  now  be  solved  for  n  =  2. 

f'(t)  =  -  rf  (t)  +  r  max  (p  (2  -  e"r  )  ,  p  ) 
2         2  12 

-rf  1       Pl 

Since  p  >  p  (2  -  e   )  when  either  t  <  -  In  or  if  p  >  2p,, 

2  V  r    2prp2         2    1' 

the  solution  to  the  differential  equation  is: 

-rt  1       Pl 

f  (t)  =  P0  (1  -  e  rc)   for  t  *  -  In or  if  p  a  2pn 

2      2  r     2p  -p         2    1 


-rt  1. 

f  (t)  =  2p,  -  e    (rPlt  +  p.  +  prt  -  p.  In 


2       1   "    v  rl    rl  '  v2       *\   -  2prp2  ' 

1       pl 
for  t  ^  -  In  and  p  <  2p„ 

2P  -p         2   Pl 

1  2 

The  significance  in  this  solution  lies  in  the  "cross-over"  point  in 
time.   If  the  search  time  remaining  is  less  than  this  value  of  t,  the 
optimal  return  is  achieved  by  dropping  both  bombs  on  the  next  target 
which  is  encountered.   If  greater,  only  one  bomb  is  dropped  on  the 
target,  the  other  being  saved  for  possible  future  encounters. 
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Having  found  f  (t) ,  f  (t)  can  now  be  calculated.  The  differential 

equation  however  becomes  unduly  complicated  for  n  >  2.  A  numerical 

solution  may  be  obtained  for  the  problem  though  by  recourse  to  equation 

(1),  the  difference  equation,  with  which  f  (t)  can  be  calculated  for 

n 

any  n  to  any  desired  degree  of  accuracy  by  taking  At  small  enough. 

With  the  possibility  of  encounters  with  targets  of  different 
types,  the  value  of  a  target  which  is  met  must  be  considered  when  the 
decision  is  made  as  to  how  many  weapons  should  be  expended. 

The  methods  used  in  developing  the  functional  equation  are 
analagou6  to  those  for  the  simpler  model  with  the  addition  of  con- 
siderations of  target  value  and  different  kill  probabilities.   Define: 

m  °  the  number  of  different  target  types. 

p    =  the  probability  that  j  weapons  expended  on  a  target 
J    type  i  will  achieve  a  kill. 

r  ■  average  encounter  rate  with  type  i  targets, 
i 

v  »  value  of  target  type  i. 

f  (t)  ■  expected  value  of  targets  killed  in  the  interval  (0,  t) 
given  n  bombs  at  time  t  =  0  and  optimal  utilization  of 
these  weapons. 


Then: 

m  m 

(3)    f  (t)  -  (1  -  )  r.  At)  f  (t  -  At)  +  '  r.  At  max  (vjp  +f   (t-At)) 

i=l  i-1      O^j^n 


Note  that  in  this  case  the  maximization  is  over  the  range  0  ^  j  :£  n 
vice  1  =£  j  ^  n.   Here  the  expected  return  might  be  greater  if  no 
weapons  are  dropped  on  targets  of  little  value,  expecting  at  some 
later  time  in  the  mission  to  meet  something  of  greater  worth. 


17 


The  differential  equation  is  obtained  by  subtraction  of  f  (t) 
from  both  sides  of  the  above  equation,  division  by  At  and  taking  the 


limit  as  At  approaches  zero. 


m 

V 


m 
V 


f,(t)  =  'I  'iVO  +  Z  ri  maX  (Vij  +  fn-j(t)) 
i=l         1-1   °^n 


For  n  -   1  and  m  =  2  the  solution  to  this  equation  can  be  shown  to  be: 
(r1P1v1  +  r2p2v2) 


f   (t) 


(ri  +  r2> 


(1   -   e"*rl  +  r2)t:  )    for  p  v  >  p  v 

2  2       1   1 


and 


(Kt^ 


1 1ti   <rlPlvl  +  r2p2v2> 

r9(P2V2   "   P1V1> 


In 


"l  +  r2 


or   if     p.v,>  p  v     then  for 
11       2  2 


A                   1            ,      <riPlVl  +  V2V 
0<st^     In =-=- 


r     +  r«  /  \ 

1         2  Wl   •  P2V 


or   if  p  v     =  p  v     then  for  all  t. 
11  2  2 


r  riPivi  +  r?Povo 

(t)  =  p  v  +  (p  v  -  p  v  )        l  l  r  +r 

11     2  2    l   X  L  r,  (P  v  -  P  v  )  J   1  2 


\  ^PlVl  "  V 


+r     e 
2 


v 


r  pnv,  +  r  p  v 
1   ,    V  1  1    2H2  2 
for  p  v  >  p  v  and  t  ^  — q: —  in  

1122  rl  r2  ri(plvi   "   P2V 


r  ripivi  +  r2P2v2 

fx(t>    =  P^  +   (p^   -P2v2)    L  —J——— 


ri+^2  2 

e 


2   XK2  2       rH' 


1  rlplvl  +  rop2v2 

for  p  v  >  p,v,    and   t  ^     — - —  In  -1 * *-* — 

2  2     "11  r  +r  r     (p  v     -  p  v.) 

12  2       2  2         1  l 
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Thus  for  the  simple  case  of  two  target  types  and  a  single  weapon, 

the  expected  returns  can  be  calculated  analytically.   In  addition,  the 

optimal  strategy  is  found,  here  consisting  merely  of  a  cross-over  time, 

the  time  before  which  it  does  not  pay  to  fire  the  weapon  at  the  target 

with  the  lowest  p  v  .   For  n  >  1  however,  it  can  readily  be  seen  that 
i  i 

solutions  to  the  differential  equation  are  not  so  easily  obtained. 
Again  though,  through  recourse  to  the  difference  equation,  (3),  numer- 
ical solutions  can  be  obtained  for  any  n  to  any  desired  degree  of 
accuracy. 

Further  generalization  of  the  model  entails  a  relaxation  of  the 
restriction  to  a  single  weapon  type.   The  functional  equation  is 
similar  to  those  for  the  simpler  models  with  some  minor  changes  in 
the  definitions  of  the  terms.   Define: 


p    ■  the  probability  that  k  weapons  of  type  j  expended  on 
*jk   a  target  of  type  i  will  achieve  a  kill. 

L  =  the  number  of  different  weapon  types. 

f  (t)  =  expected  value  of  targets  killed  in  the  Interval 

n,n2n3",nL      (0,  t)  ,  given  n.  weapons  of  type  i  (for  1*1,  •••,L) 
at  t  ■  0  and  optimal  utilization  of  these  weapons. 

The  difference  equation  can  be  written: 

in 

f  (t)  =  (1  -V  r  At)  f  (t  -  At) 

n  n  n  ...n  u     i     n,non-,"*nT 

1  2  3    L        i=i        1  2  3    L 


m 

+  )    r  At     max         max  (v.p         +  f  (t   -  At)) 

L*     i  x   ijk         n  n   ...n  -k...n 

i=1  l*j£L     0<fc<n  1  2  j  L 
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The  differential  equation  can  be  written  as  before,  but  the 
solutions  are  obtainable  only  in  the  most  trivial  cases. 

A  FORTRAN  program  has  been  written  however  from  which  numerical 

solutions  to  the  difference  equation  can  be  calculated.   With  inputs 

of  the  aircraft  load  carrying  capabilities,  the  mission  time,  and  the 

parameters  of  the  model  (m,  r.,  v  ,  L,  p..,  ,  and  At)  the  program 

i   i      ijk 

selects  values  for  the  n  which  are  consistent  with  the  limitations 

i 

of  the  aircraft  and  maximize  f .  Additionally,  another  FORTRAN  pro- 
gram lists  the  optimal  strategy  to  be  employed  in  the  decisions 
relating  to  the  expenditure  of  these  weapons. 


A  SUB -OPTIMAL  DETERMINATION  OF  LOAD 

In  the  majority  of  practical  cases  hand  calculation  will  be  the 
only  means  available  for  determining  loading  requirements.   In  such 
cases  a  sub-optimal  method  can  be  used  at  the  squadron  level  for 
obtaining  the  initial  weapons  load  for  a  mission.  Define: 


T  =  search  time  (i.e.,  mission  time  allotted  to  searching  for 
and  attacking  targets). 

P   (T)  =  probability  that  target  type  i  will  be  met  n  times  during 
the  mission. 

00 

n.  (T)  "  /  P   (T)  *  probability  that  target  type  i  will  be  met  at 
n=^  ni      least  once  during  the  mission. 


v  =  value  of  target  type  i. 
i 

m  =  number  of  target  types. 

b   =  number  of  sticks  of  s  weapons  of  type  j  in  the  load, 
js 

L  ■  number  of  weapon  types  available. 


■•■In  this  development  the  terms  stick  and  salvo  are  used  equiva- 
lently.  The  choice  of  a  stick  or  a  salvo,  in  the  load  determination, 
will  depend  upon  which  gives  the  higher  kill  probability. 
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k    ■  probability  that  a  stick  of  s  type  w  weapons  will  kill 
iws   target  type  i. 

Note  that  in  this  development  individual  bombs  are  considered  as 
a  stick  of  one;  i.e.,  two  250  lb.  bombs  might  be  considered  as  two 
sticks  of  one  type  j  or  one  stick  of  two  type  j  weapons. 

If,  a  priori,  only  the  more  highly  valued  targets  and  those 
individual  weapons  with  the  highest  kill  probabilities  against  these 
targets  are  considered,  the  assumption  can  be  made  that  a  target 
sighted  will  be  dropped  on,  as  long  as  the  appropriate  weapon  for 
that  target  is  available.   (Note  that  this  assumption  is  used  only 
for  determining  the  load.   In  the  actual  performance  of  the  mission 
there  is  the  possibility  that,  if  that  weapon  most  effective  against 
a  target  is  not  available,  a  less  effective  weapon  will  be  dropped.) 
The  problem  then  becomes  one  of  maximizing,  for  a  given  mission  time, 
the  expected  value  of  targets  killed,  E,  for  those  loads  being  con- 
sidered.  E  may  be  determined  as  follows: 

m 

j.Vf    k    v 
E  L     ijs   ijs   i 

i=l 

where  f    is  the  expected  number  of  type  i  targets  on  which  weapon 
ijs 

type  j  can  be  expended,  and  j  indicates  that  weapon  type  for  which 

k, u  *  k    for  all  w  =  1,  2,  *•*,  L. 
■LJX    iwl 

If  the  number  of  type  j  weapons  available  were  limitless,  f  . 

oo  ijs 

would  simply  be  equal  to  ri  =  ,  np  .  (T)  .   However,  since  a  given  load 

n=o 

is  limited  to  b   sticks  of  type  j  weapon,  no  more  than  b.   can  be 
js  js 

dropped  regardless  of  how  many  type  i  targets  are  sighted.   Thus 


The  decision  as  to  which  targets  and  weapons  will  be  considered 
will  depend  highly  upon  the  judgement  of  the  individual  responsible  for 
determining  the  aircraft  load. 

21 


fijg  is  given  by 


b.  -1 
js 

(4)       £,    -^  up  .(T)  +  £   b  p  .  (T) 

ijs   /-■   ni      *-'    js  ni 

n=o        n=b 

js 

Equation  (1)    can  be  reduced   to  a  more  easily  calculable   form. 


SinCe  bjs-l  bjs-l  co  b      -1 

iVd)    =1     ^(T)    and     V       Pn.(T>    =  H  .-     J     p^fc) , 

n=o  n=l  n=b.  n~-L 

js 

equation   (4)    can  be  written  as 

bjs-l  bjs-l 

£     -;    np    (t)  +b.  n  -  b.     I  pni(T> 

1JS     ^         ni  js   i  js  , 

n=l  n-x 


b.    -1 
js 

=  b     n     -Y  (b.  -n)  p     (I) 

js  i      ^     Js  ni 

but  II.    =   1  -   p      (T)    and  consequently     b     II    =  b     -   b.   p      (T) .      Note 
1  Oi  js   i  js        js  Oi 


that  b     p        can  be  expressed  as   (b     -  0)    p      .     Thus  equation  (4) 
js  Oi  Js  Oi 

now  reduces   to  ,         , 

b  - 1 

js 

(5)       f        =  b   -  )       (b -n)  p   (T) 
ijs   js   ^    js     ni 
n=o 

Equations  (4)  and  (5)  are  valid  for  all  b   ^  1.  By  definition 

js 

f .  .  =0  for  b.   =0,  since  no  sticks  of  type  j  weapon  can  be  dropped 
if  none  are  carried. 

As  noted  before,  the  random  distribution  of  targets  and  conse- 
quent exponential  inter- encounter  times  lead  to  a  Poisson  counting 
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process  in  determining  the  probabilities  of  encountering  given  numbers 

of  targets;  i.e.,  for  type  i  targets  appearing  at  a  mean  rate  r  , 

i 

p  .(T)  -  e"riT  £#n  . 
ni  n! 

This  gives  D   -1 

JS  T  /  Tin 

-r  T  (riT) 

(6)       f    =  b   -    (b1s-n)  e  i  ^ 

ijs    Js   <-    Js  n- 

n=o 

Table  11,  Appendix  III  gives  values  of  f..   for  various  com- 

i-JS 

binations  of  b   and  r  T.   This  table  greatly  facilitates  calculation 
of  a  "best  load"  from  the  feasible  loads  being  considered. 

To  determine  the  "best  load"  for  a  given  set  of  targets,  find 
that  individual  weapon  type  j,  for  each  target,  which  has  the  highest 
kill  probability.   The  loads  considered  are  then  the  different  mixes 
of  these  weapons  selected  from  a  list  of  feasible  loads  for  the  air- 
craft. 

For  a  given  load,  b.   will  be  known  for  each  type  weapon  in  the 

J  s 

load  and  the  encounter  rate  r  will  be  estimated  for  each  of  the 

i 

critical  targets.   The  planned  search  time  T  will  also  be  known  for 

the  mission.   Thus  using  known  values  of  k. .   and  estimated  values 

1JS 

of  v  ,  E  can  be  calculated. 

Consider  first  the  case  where  weapons  will  be  dropped  indivi- 
dually (sticks  of  one  only)  and  where  each  weapon  type  in  the  load 
will  be  used  against  only  one  type  target.  In  other  words,  target 
types  and  weapon  types  are  matched  one  to  one.   Since  weapons  are  to 

be  dropped  singly,  b   e  b  ..  for  each  weapon  type  j.  To  determine  E 

js    J  l 

for  such  a  load,  first  determine  b   and  r.T  for  each  weapon- target 

J1      L 
type  combination.  Then  obtain  f    from  Table  11,  and  multiply  this 

1J6 

23 


value  by  k   v  .   Let  f   k   v  =  E     Add  the  resulting  E   's  to 
ijs  i        ijs  ijs  i    ij  ij 

m 
set  E  for  the  load,  i.e.,  E  =  ^v        .   (Note  that  the  E..'s  are  summed 

i=l 
over  i  only,  since  the  subscript  j  simply  indicates  that  weapon  which 

is  most  effective  against  target  type  i.)  In  the  same  manner,  E  may 
be  determined  for  each  load  being  considered;  that  load  which  yields 
the  maximum  E  is  then  the  "best  load." 

For  this  sub-optimal  model,  the  extension  of  the  problem  to 
include  sticks  of  two  or  more  weapons  is  accomplished  quite  easily, 
since  different  size  salvos  or  sticks  of  a  given  type  weapon  con- 
stitute different  loads.  For  example,  as  noted  previously,  two  250 
lb.  bombs  might  be  considered  as  a  load  consisting  of  two  sticks  of 
one  bomb  or  as  another  load  consisting  of  one  stick  of  two  bombs. 
With  this  in  mind,  E  can  be  calculated  for  these  loads  in  the  same 
manner  as  previously  outlined.   In  order  to  keep  the  problem  hand 
workable,  this  approach  assumes  that  each  delivery  of  a  given  type 
weapon  will  consist  of  the  same  size  stick.  This  assumption  is  made 
solely  for  determination  of  the  aircraft  loading  prior  to  flying  the 
mission.   During  the  actual  mission,  however,  it  is  quite  likely 
that  sticks  of  different  sizes  will  be  expended  at  different  times. 
For  example,  a  pilot  might  drop  one  type  j  bomb  on  the  first  truck 
he  sees,  a  stick  of  three  type  j  bombs  on  the  second  truck  he  sees, 
etc.   Prior  to  the  mission,  however,  this  load  of  four  type  i  weapons 
could  have  been  evaluated  as  one  stick  of  four,  two  sticks  of  two,  or 
possibly  as  four  sticks  of  one,  since  calculations  for  sticks  of 
varying  sizes  quickly  exceed  the  bounds  of  practical  hand  workability. 
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In  addition  to  the  case  where  targets  and  weapons  are  matched  one 
to  one,  consideration  must  be  given  to  the  case  where  a  given  weapon 
type  i6  optimal  for  more  than  one  target  type.   If  the  assumption  that 
a  target  sighted  will  be  dropped  on  were  adhered  to,  the  expected 
number  of  type  i  targets  killed  would  be  proportional  to  the  encounter 
rate  r,.   However,  in  this  situation,  a  concession  is  made  to  the  fact 
that,  early  in  the  mission,  a  target  of  lesser  value  might  be  passed 
in  hopes  of  encountering  a  more  highly  valued  target  later.   Thus  it 
is  assumed  that  the  expected  number  of  type  i  targets  killed  is  pro- 
portional to  both  the  encounter  rate,  r  ,  and  the  target  value,  v  . 

i  i 

The  expected  return  from  a  given  load  is  now  calculated  exactly  as 
before  with  the  following  minor  modification:   For  all  targets  type  i 
against  which  weapon  type  j  is  most  effective,  calculate  the  E-.'s 
exactly  as  in  the  simpler  case  (as  if  all  type  j  weapons  were  to  be 
expended  on  each  type  target).   Then  to  approximate  the  particular 
expected  value  of  type  y  targets  killed  by  weapon  type  j,  E   is  mul- 

yj 

rv  vv 
tiplied   by   the  normalizing   factor  N     =  — * — *—     where   the   summation 

Srivi 

in  the  denominator  is  taken  over  those  type  i  targets  for  which  weapon 

type  j  is  optimal.  These  expected  values  are  then  added  to  give  E  as 
before. 

A  logical  extension  to  this  hand  approximation  is  the  consider- 
ation of  "secondary"  weapons.  A  secondary  weapon  would  be  considered 
to  be  a  weapon  which  has  high  kill  probabilities  against  several 
targets,  yet  is  optimal  against  none  of  them.  There  are  certainly 
cases  where  such  a  weapon  might  give  higher  expected  returns  over  the 
entire  target  spectrum  than  would  a  load  consisting  of  only  those 
weapons  which  are  optimal  for  the  individual  targets.  For  example, 
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consider  a  case  in  which  there  are  two  targets  and  three  weapons  to 
be  used  against  these  targets.  The  single  stick  kill  probabilities 
for  the  individual  weapons  are  as  follows: 


WEAPON 

TARGET 
.09 

1 

TARGET  2 

1 

.01 

2 

.01 

.09 

3 

.08 

.08 

If  the  load  were  limited  to  those  weapons  which  are  optimal  for 
individual  targets,  only  weapons  1  and  2  would  be  considered.  However, 
it  is  obvious  that  weapon  3  merits  examination,  and,  in  fact,  it  appears 
as  if  a  pure  load  of  weapon  type  3  might  possibly  be  the  best  load  to 
use  in  this  particular  case.   Calculation  of  E  for  such  a  secondary 
weapon  is  identical  to  that  in  the  foregoing  case  in  which  one  weapon 
is  considered  for  use  against  more  than  one  target. 

This  sub-optimal  development  may,  at  first  glance,  appear  to  be 
long  and  cumbersome.   It  must  be  realized  though  that  in  most  actual 
situations  all  but  a  very  few  feasible  loads  will  be  eliminated  from 
consideration  prior  to  any  calculation,  and  these  will  be  examined 
only  as  far  as  they  pertain  to  a  small  number  of  critical  target  types 
(normally  not  more  than  three  or  four).  Thus  it  can  be  seen  that 
relatively  few  calculations  will  actually  be  required  and  these  will 
be  of  a  trivial  nature. 

In  order  to  clarify  the  use  of  the  sub-optimal  model,  the  follow- 
ing simple  example  is  presented  in  detail. 

Consider  an  armed  reconnaissance  mission  with  a  planned  search 
time  of  one  hour.   (Thus  T  =  1  so  that  r.T  =  r.  for  all  targets). 
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Three  target  types  are  considered: 


v  «=  2 
1 


v  -  1 

2 


v  =  4 
3 


ri  =  5 


r  =  3 
2 


r3  =  2 


The  aircraft  has  two  weapon  stations  and  two  types  of  weapons 
are  considered,  one  more  effective  against  target  type  1,  the  other 
more  effective  against  types  2  and  3  (see  Table  2). 

The  possible  loads  are  indicated  in  Table  1. 


! 

WEAPON 

1 

2 

Load 

Single 

Stick  of  2     Single 

Stick  of  2 

A 

1 

1 

B 

2 

C 

2 

D 

1 

E 

1 

TABLE  1 
COMPOSITION  OF  POSSIBLE  LOADS 


s   i 


.02   .03   .01 


.06   .08   .03 

I 1 


.01   .05   .02 


.03   .12   .05 


TABLE  2 
KILL  PROBABILITIES 
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Load  A:   b   =  1,  b   =1 
11       21 

Target  1  -  Weapon  1 

f    =  .99326  (from  Table  11;  b .  =  1,  r  T  =  5) 
111  Js      i 

E    =  f   k   v  =  (.99326)  (.02)  (2)  =  .0397 
11     111  111  1 

Target  2  -  Weapon  2 

f    =  .95021  (b   =  1,  r  T  =  3) 
221  js       i 

E    =  f   k   v  =  (.95021)  (.05)  (1)  -  .0475 
22     221  221  2 

Target  3  -  Weapon  2 

f    =  .86466  (b   -  1,  r  T  -  2 
321  js       i 

E    =  f   k™-,v  =  (.86466)  (.02)  (4)  =  .0692 
32     321  321  3 

Since  weapon  type  2  is  optimal  for  two  target  types,  E   and  E 

must  be  normalized. 

N  E       =   (r  v  /r  v     +  r  v  )    E       =  (3/11)    (.0475)    -   .0129 
2  22  2222         33       22 


NE       =  (8/11)    (.0692)   =   .0608 
3  32 

Thus  for  load  A,  E  =  .0397  +  .0129  +  .0608  -  .1134 

Load  B:   b   =2 
11 

f    =  1.95283  (b.  =  2,  r  T  =  5) 
111  js       i 

E    =  E  (Load  B)  =  (1.95283)  (.02)  (2)  =  .0724 

Load  C:   b   =2 
21 

f    =  1.75106  (b.  =  2,  r  T  =  3) 
221  J       i 
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N2E21  "  (3/U)  (1"75106>  (-05)  <L)  =  -02^° 

f     »  1.45866  (b   =  2,  r  T  =  2) 
321  Js       L 

N  E   =  (8/11)  (1.45866)  (.02)  (4)  =  .0846 
3  32 

E(Load  C)  -  .0240  +  .0846  =  .1086 

Load  D:   b   =1  (Note  that  this  is  almost  the  same  as  Load  B,  but 
12 

now  the  two  weapons  are  considered  as  one  stick  of  two  vice  two  sticks 
of  one) . 

f    =  .99326 
112 

E  -  E(Load  D)  -  (.99326)  (.06)  (2)  =  .1192 

Load  E:   b   =1 
22 

f    =  .95021 
222 

N  E   -  (3/11)  (.95021)  (.12)  (1)  -  .0312 

2  22 

f    =  .86466 
322 

N  E   =  (8/11)  (.86466)  (.05)  (4)  =.1252 

3  22 

E  (Load  E)  =  .0312  +  .1252  =  .1564 

For  the  loads  considered,  the  maximum  E  is  .1564.  Therefore  the 
"best  load"  is  load  E  which  consists  of  one  stick  of  two  type  2  weapons 
Note  that  if  sticks  of  one  only  (case  1  in  the  model  development)  had 
been  considered,  the  "best  load"  would  have  been  load  A  with  one  of 
each  type  weapon  yielding  E  =  .1134. 

The  foregoing  sub-optimal  method  is  limited  in  that  it  considers 
search  time  as  a  whole  rather  than  a  time-stepped  firing  doctrine  as 
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in  the  dynamic  programming  method.  Thus  those  added  values  which 
might  be  generated  by  strategies  which  vary  with  mission  time  are 
lost  in  this  formulation.   Consequently  the  expected  kill  values 
will  generally  be  much  lower  than  those  of  the  dynamic  program  formu- 
lation, and  the  "best  load"  generated  will  not  always  be  the  optimum 
load.   However,  the  load  obtained  by  this  method  should,  if  expended 
in  accordance  with  the  optimum  firing  doctrine  of  the  dynamic  program, 
be  one  of  the  more  highly  valued  loads. 
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3.   SAMPLE  RESULTS 

OPTIMAL  LOAD  AND  FIRING  STRATEGY  DETERMINATION 

For  a  mission  time  of  15  minutes,  seven  weapons  from  which  to 
choose,  and  values  of  the  other  input  parameters  as  listed  in  Table 
3,  an  optimal  loading  has  been  determined  using  the  "BESTLD"  program 
of  Appendix  II,  with  the  following  inputs: 

Weapon  Type  :    1     2     3     4     5     6      7 

Weapon  Weight:      265        300       510        600        675        950        1975 

Number  of   racks   ■  5 

Max  pay load   -  5,000    lbs. 

Number  of   target   types   ■  4 

Mission   time   -  0.25   hours 

Delta  T   ■    .010  hours    (This    is    the  At  of   equation   (3)) 


TARGET 

TYPE 

ENCOUNTER  RATE 
30.00000 

VALUE 
1  .00 

\ 

1 .00000 

2  5.00 

5 

.  10003 

6  0.00 

3.  Of 

5.00 

KILL 

PROBABILITIES 

NUM 

OF  WE 

=     1 

2 

3 

4 

5 

TARGET 

1 

WEAPHN 

1 

.02C 

.040 

.06  0 

.070 

.080 

TARGET 

i 

WEAPON 

2 

.050 

.QUO 

.06  0 

.070 

.0  80 

TARGET 

1 

WEAPON 

3 

.0  40 

.0  70 

.  09  0 

.  100 

.1  10 

TARGET 

1 

..  'PON 

4 

.00  5 

.00  8 

.010 

.011 

.012 

TARGET 

i 

WEAPON 

5 

.  0  3  C 

. 

.(-70 

.071 

.0  80 

TARGET 

i 

WEAPON 

6 

.05  0 

.0  80 

.10  0 

.  120 

.130 

TARGET 

1 

W EAPON 

7 

.  120 

.  160 

.  180 

.200 

.202 

TARGET 

t 

WEAPON. 

1 

.  0  1  0 

.0  20 

.02  4 

.030 

.0  35 

TARGET 

WEAPON 

2 

.020 

.030 

.04  0 

.050 

.060 

TARGET 

2 

WEAPON 

3 

.012 

.021 

.  03  0 

.0  36 

.0  42 

TARGET 

2 

WEAPON 

4 

.016 

.0  32 

.04  7 

.061 

.076 

TARGET 

-■ 

WEAPON 

5 

.02  2 

.039 

.053 

.065 

.07  5 

TARGET 

•i 

WEAPON 

6 

.007 

.012 

.  0 1  6 

.019 

.02  1 

TARGET 

2 

WEAPON 

7 

.  OS  G 

.  150 

.  153 

.200 

.240 

TARGET 

j 

WEAPON 

1 

.0  50 

.080 

.  140 

.170 

.190 

TARGET 

-. 

WEAPON 

2 

.094 

.  1U3 

.172 

.  196 

.20  5 

TARGET 

•j 

WEAPON 

3 

.06C 

.090 

.  150 

.190 

.22  0 

TARGET 

b 

wEAPOM 
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.036 

.0  70 

.102 

.131 

.161 

TARGET 

3 

WEAPON 

5 
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.  175 

.22  0 

.233 

.245 

TARGET 

3 

WEAPON 

6 

.07C 

.  100 

.  160 

.200 

.230 

TARGET 

3 

WEAPON 

7 

.08  0 

.  120 

.180 

.220 

.250 

TARGET 

n 

WEAPON 

1 

.01  1 

.020 

.028 

.035 

.042 

TARGET 

4 

WEAPON 

2 

.019 

.036 

.04  8 

.056 

.064 

TARGET 

k 

WEAPON 

3 

.019 

.032 

.042 

.05  1 

.057 

TARGET 

h 

WEAPON 

U 

.00  1 

.001 

.001 

.001 

.002 

TARGET 

k 

WEAPON 

5 

.052 

.054 

.05  4 

.082 

.090 

TARGET 

k 

WEAPON 

6 

.00  9 

.016 

.02  3 

.028 

.033 

TARGET 

4 

WEAPON 

7 

.029 

.057 

.066 

.080 

.090 

TABLE   3 
KILL  PROBABILITIES 
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LOAD  NO. 

LOADING 

1 

0300002 

2 

1200002 

3 

0100102 

4 

1000102 

5 

3000002 

381 

2011010 

These  weapons'  weights  and  aircraft  payload  and  rack  capacities 
result  in  approximately  400  possible  loadings  which  must  be  considered. 
From  amongst  these  400,  Table  4  lists  the  five  loads  with  the  greatest 
expected  returns  and  for  comparison  the  loading  yielding  the  least. 
A  loading  designated  1200002  denotes  one  weapon  of  type  one,  two 
weapons  of  type  two,  zero  weapons  of  types  three  through  six,  and  two 
weapons  of  type  seven. 

EXPECTED  RETURN 

1.29514 

1.25048 

1.22353 

1.20670 

1.20409 

.38069 

TABLE  4 
LOADINGS  AND  EXPECTED  RETURNS, 
DYNAMIC  PROGRAMMING  METHOD 

The  expected  returns  from  the  optimal  load  are  over  three  times 
higher  than  those  of  the  worst.   In  fact  the  twentieth  best  load, 
0020300,  had  expected  returns  of  only  .77700.  Thus  the  optimal  loading 
yields  returns  more  than  65%  better  than  a  loading  which  is  better  than 
9570  of  the  possible  loads. 

The  above  inputs  and  constraints  have  also  been  used  in  the  sub- 
optimal  model  to  generate  a  "best  load."  Weapon  7  is  the  best  single 
weapon  against  targets  1  and  2,  and  weapon  5  is  best  against  targets 
3  and  4.  Thus  the  only  loads  considered  in  this  formulation  were 
0000500,  0000401,  and  0000102.  Table  5  lists  six  possible  load  con- 
siderations (including  stick  options)  and  the  expected  returns  from  each. 
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Loading  E 

0000102  (Weapon  7  used  as  a  stick  of  2)  .5766 

0000401  (Weapon  5  used  as  a  stick  of  4)  .4458 

0000102  (Weapon  7  used  as  2  sticks  of  1)  .3824 

0000500  (Weapon  5  used  as  a  stick  of  5)  .3740 

0000401  (Weapon  5  used  as  4  sticks  of  1)  .3320 

0000500  (Weapon  5  used  as  5  sticks  of  1)  .1432 

TABLE  5 
POSSIBLE  LOADS  AND  EXPECTED 
RETURNS,  SUB-OPTIMAL  METHOD 

The  "best  load"  using  the  sub-optimal  hand  method  is  0000102 
considering  weapon  7  as  one  stick  of  two.  Although  only  the  principal 
weapon  for  each  target  lias  been  considered,  note  that  weight  constraint 
still  allows  the  loading  of  one  more  of  the  lighter,  less  effective 
weapons  (weapon  2  rather  than  weapon  1  since  weapon  2  is  dominant  over 
the  entire  target  spectrum) .   The  load  thus  generated  corresponds  to 
the  third  most  valuable  dynamic  programming  load,  and,  if  used  in 
accordance  with  the  dynamic  programming  firing  policy,  it  would  be 
expected  to  yield  a  return  of  greater  than  907,  of  that  given  by  the 
actual  optimum  load. 

With  the  best  loading  as  generated  by  either  method,  the  optimal 
strategy  can  be  listed  using  the  "POLICY"  program  of  Appendix  II.  The 
results  of  this  determination  for  loading  0300002  are  shown  in  Figure 
1  along  with  a  format  for  displaying  the  strategy  compactly.  The 
time  axis  is  mission  time  remaining  and  entries  such  2x7  denote  that 
two  weapons  of  type  seven  should  be  expended. 

As  an  example,  as  the  mission  begins,  the  subloading  is  0300002. 
If  the  first  target  which  is  met  is  a  type  one,  and  is  encountered  with 
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FIGURE  1 
OPTIMAL  STRATEGY  DISPLAY 
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(The  abscissa  in  these  figures  is  the  time  remaining  in  which  to 
search  for  and  attack  targets.) 


36 


.22  hours  remaining,  the  optimal  strategy  indicates  that  no  weapons 
should  be  fired.   If  the  next  target  which  is  encountered  is  a  type 
four  at  time  .17,  the  subloading  still  being  0300002,  a  drop  of  three 
weapons  of  type  two  is  called  for.   The  sub-loading  now  is  0000002. 
If  a  type  three  target  is  now  met,  both  of  the  remaining  type  seven 
weapons  should  be  dropped. 

It  is  of  interest  to  note  the  changes  which  occur  over  time  in 
the  optimal  weapon  to  drop  when  a  target  is  met.   With  the  loading 
0300002  and  target  type  one,  with  more  than  .08  hours  remaining,  no 
weapons  should  be  expended.   With  .045  to  .08  hours  remaining,  the 
optimal  strategy  calls  for  tne  expenditure  of  one  weapon  of  type  two. 
Weapon  type  seven  however  lias  a  higher  kill  probability  against  type 
one  targets'   The  type  seven  weapons  at  this  point  are  being  saved 
for  possible  future  encounters  with  targets  of  type  two  and  four,  botn 
of  which  are  more  valuable  and  occur  frequently  enough  so  that  there 
is  a  significant  probability  of  an  encounter  in  the  time  remaining. 
Although  weapon  ty^e  two  is  the  most  potent  against  type  three  targets 
-the  most  valuable  target  type-  the  probability  of  meeting  a  type 
three  target  is  rather  small  because  of  the  low  encounter  rate  of  0.10, 
It  is  for  this  reason  that  type  two  weapons  are  not  being  conserved. 
Later  in  the  mission  however,  the  optimal  weapons  to  drop  on  a  type 
one  target  are  type  sevens.  At  this  point  in  time  the  probability  of 
future  encounters  with  more  valuable  targets  is  remote. 

Since  the  ultimate  objective  is  the  maximization  of  the  expected 
returns  for  the  mission,  a  consideration  might  be  made  for  lengthening 
the  search  time  at  the  expense  of  lower  payload  and  weapon  station 
availability.   That  is,  perhaps  an  additional  tank  of  fuel  can  be 
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carried  on  one  of  the  aircraft's  weapon  racks,  thus  reducing  the 
number  available  on  which  ordnance  can  be  carried  and  also  reducing 
the  weapons  payload,  but  increasing  the  search  time. 

Table  6  lists  the  expected  returns  for  the  five  best  loads 
generated  when  the  search  time  has  been  increased  to  0.50  hours,  the 
payload  decreased  to  4000  lbs.,  the  number  of  weapon  racks  decreased 
to  four,  and  the  other  input  parameters  unchanged  from  those  of  Table 
3. 

LOADING     EXPECT ED  RETURNS 

0000002  1.78303 
0000301  1.42184 
0100201  1.30025 
0010231  1.29748 
0000400  1.28808 
TABLE  6 
BEST  LOADS  AND  EXPECTED  RETURNS 

Since  the  expected  value  of  targets  killed  for  the  longer  mission 
is  significantly  better  than  that  for  the  shorter,  1.78803  compared  to 
1.29514,  it  is  clearly  preferrable  to  make  the  trade-off  of  gas  for 
weapons. 

An  interesting  feature  of  the  kill  probability  matrix  of  Table 
3  is  the  preferred  weapons  with  which  to  attack  the  most  valuable 
targets.  Weapon  type  seven  is  the  most  potent  against  type  two  tar- 
gets while  weapon  type  five  is  preferred  against  type  three  targets. 
Target  type  three,  although  the  most  valuable,  has  a  very  low  encounter 
rate.  This  results  in  the  fact  that  weapon  type  five  does  not  appear 
frequently  in  the  loadings  which  yield  high  expected  returns  for  short 
missions.   Letting  the  maximum  payload  equal  4000  lbs.  and  the  number 
of  weapons  stations  equal  four,  the  preferred  loadings  of  these  two 
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weapons  alone  can  be  examined  as  a  function  of  search  time.  The  three 
loadings  which  must  be  considered  are  02,  31,  and  40,  where  now,  02 
denotes  a  load  consisting  of  zero  weapons  of  type  five  and  two  weapons 
of  type  seven.   Calculations  show  that  an  initial  loading  of  02  is 
optimal  for  mission  times  of  less  than  1.6  hours.  For  mission  times 
of  1.6  to  8.9  hours,  loading  31  is  optimal,  and  finally  for  T  >   8.9 
hours,  loading  40  is  preferred. 

COMPARISON  OF  EXPECTED  RETURNS  FROM  NON-OPTIMAL 
FIRING  STRATEGIES 

The  gains  in  kill  potential  realizable  through  the  utilization 

of  the  optimal  firing  doctrine  can  be  demonstrated  by  comparisons  with 

the  expected  returns  from  other  strategies.   As  an  example,  consider 

the  following  case  with  non-optimal  strategies: 

Let: 

Number  of  weapons  =  1 

Operating  time  ■  1  hour 

No.  of  target  types  =  2 

v  ■  3 

v  =4 
2 

Strategy  A:   Expend  the  weapon  on  the  first  target  which  is 
encountered. 

E.  =  the  expected  returns  from  strategy  A. 
A 

Strategy  B:   Expend  weapon  only  if  a  target  of  type  2  is 

encountered.   (Note  that  target  type  2  is  the 
more  valuable)  . 

E  =  the  expected  returns  from  strategy  B. 
B 

T.:   The  random  variable  describing  the  time  to  first  encounter 
with  target  type  i.   (T .  is  exponentially  distributed  with 


mean  1/r . .) 

i 


L 


39 


E  =  p  v..Prob.  (T  <  min  (t  ,  1  hour))  +  p  v  Prob.  (T  <  min  (t  ,  1  hour)) 
A*1        l  2  22        2        1 


E  =  PnV 


-r  t  °°,    -r  t 
re  l  l       re  2  2  dt  dt  +  p0v 
1      ^   2        2  1    2  2  J 
t 

1 


r  e 
2 


-r  t 
2  2 


-r  t 

re  1  1  dt  dt 

*  l  l     2 

t2  * 


E  -  p  v  Prob.  (T  <  1  hour) 
B    2  2        2 


-r  t 
E  ■  P^v^   r  e  2  2  dt 


B 


2  2  J 


For  several  values  of  the  parameters  p.  and  r  ,  Table  7  shows  the 

1      i 

results  of  calculations  of  the  expected  returns,  using  either  of  these 
strategies  or  the  optimal. 


Expected  Value 

of  Targets  Killed 

%  Gain  Optimal 

Strat. 

Strat. 

Optimal 

Over 

A 

B 

Strat. 

A                  B 

v-3 

P2  *  -4 

r     -   10 
1 

.963      i  1.011 

1.270 

32% 

26% 

r     =  1 

2 

\ 

Px-  .3 

!                        1 

1 

1  r     =  10 
1 

! 
.916            .354 

1.018 

11% 

185% 

r     =   .25 

'  2 

If    -  .1 

• 

i 

1 

P2  =  .4 

.414         1.011       1.064 

150% 

5% 

rl  "  10 

i 

!  r     =  1 
2 

. > 

! 

TABLE  7 
GAINS  IN  EXPECTED  VALUE  OF  TARGETS  KILLED  ACHIEVED  BY 
UTILIZING  THE  OPTIMAL  FIRING  STRATEGY 
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These  results  are  representative  of  those  obtainable  using  other 
vciiues  for  the  input  parameters.   Generally,  when  one  of  these  non- 
optimal  strategies  approaches  the  optimal  in  expected  returns,  the 
other  yields  very  poor  results.   When  the  non-optimal  strategies  yield 
roughly  equivalent  returns,  the  optimal  dominates  both  by  a  significant 
amount;  this  is  demonstrated  in  the  first  example  of  the  table  by  gains 
of  26%  and  32%  through  use  of  the  optimal  strategy. 

Comparisons  can  be  made  in  this  way  with  any  conceivable  strategy. 
It  seems  reasonable  too  that  techniques  which  consider  the  target 
values  and  encounter  rates  could  yield  returns  comparable  to  those 
achievable  through  use  of  the  optimal  strategy.   The  determination  of 
these  strategies  however  would  probably  involve  calculations  no  less 
complex  than  those  of  the  dynamic  program. 

SENSITIVITY  OF  EXPECTED  RETURNS  TO  TARGET  VALUES 

A  complete  parametric  analysis  is  beyond  the  scope  of  this  paper; 
however,  some  further  comment  must  be  made  about  target  values.   The 
developments  of  both  the  dynamic  programming  model  and  the  hand  appro- 
ximation are  based  upon  the  assumption  that  a  ratio  scale  of  values  can 
be  assigned  to  all  possible  targets.   Such  a  determination  is  completely 
subjective  and  must  rely  solely  on  the  judgement  of  experienced  planners. 
Quantitative  target  values  might  be  determined  by  presenting  a  list  of 
targets  to  several  military  commanders  and  asking  them  to  rate  these 
targets  according  to  military  value.   Targets  could  be  assigned  numbers 
consistent  with  their  importance  relative  to  a  target  of  least   value. 
The  results  of  such  a  large  number  of  opinions  could  then  be  analyzed 
by  methods  of  sensory  psychophysics  in  order  to  assign  representative 
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values  to  all  targets.  Methods  for  accomplishing  such  an  analysis 
have  been  reviewed  and  consolidated  by  S.  S.  Stevens. 

The  above  method  will  obviously  generate  strategic  values.  These 
will  certainly  not  be  valid  for  every  tactical  situation  (e.g.,  a  tank 
about  to  overrun  friendly  troops  is  more  important  than  a  railroad 
bridge  at  that  moment  even  though  the  bridge  may  have  a  higher  strategic 
value).  A  strategic  scale  will,  however,  provide  a  starting  point  -- 
a  base  which  can  be  adjusted  to  meet  specific  tactical  situations. 

An  analysis  of  the  effects  on  the  expected  returns  of  incorrect 
selections  of  target  values  can  be  made  rather  easily  for  the  simple 
case  of  a  single  weapon  and  two  target  types.  Assume  the  "true"  ratio 

of  the  target  values  is  V  *.   Now  suppose  that  an  erroneous  ratio, 

r 

V  ,  has  been  used  to  calculate  the  strategy.  The  expected  returns 

which  result  from  the  utilization  of  this  strategy  -  these  returns 

being  computed  using  the  true  value  ratio,  V  *  -  can  be  compared  with 

r 

the  maximum  of  the  expected  returns,  this  maximum  being  obtained  when 

the  strategy  derived  using  V  *  is  utilized.  Define: 

r 

vr  ■  y  vi 

and   let   the  mission  time  =   1  hour.     With  a   single  weapon,    a   strategy 
consists  merely  of  a   cross-over  time,    the   time  after  which  the   least 
valuable   target   should  be  attacked   if  met.     These  times  can  be  calcu- 
lated from: 

rlPl  +  r2P2Vr 
t L_  ln     LL    .  l  r      if  V     >  1 


Vr2         r2<P2Vr  -Pi> 


rlPl  +  r  p2V 
t  »  — L_  In  — *_£_£     if  v     <  1 

rl+r2         ri(Pl  -  p2Vr) 

t  «  "  if  V     =   1 

r 


•'-S.S.  Stevens,  "A  Metric  for  Social  Concensus,"  Science,  vol.  151, 
4  February  1966,  pp.  530-541. 
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Using  these  cross-over  times  then,  the  expected  returns,  E,  are  cal- 
culated from: 

%        "<rl+r  2>     Plrl        "<rl4*2) 

E  '  PoV  — IT"  <X  "  e       >  +  777"  (1  "  e      Z  )  if  t  U  hour 
2  r  r  +r  r  +  r 

12  12 

-r  (1-t)    -r2(l-t)      -(r  +r  )t   r 

E-pV*(l-eZ      +  e        (1-e        )   —±—   ) 

2  r  Vr2 

-r  (1-t)      -(r.+r)t   ri 
+  P  e   2     (I   -   e        l     2     \   — —     if  t  <  1  hour,  and  V  >  1 
Hle        K  e        '  r  +r  r 

1  2 

-r,(l-t)    -r^Cl-t)      -(r  fr2)t   r, 
E  -  p,(l  -  e       +  e        (1-e 


1  rl+r2 


-r  (1-t)       -<r  +r  )t   r 

+  p  V  *  e         (1-e      Z)  — - —  if  t  <  1  hour,  and  V  <  1 
P2  r  v  '  r  +r2  r 


Letting: 

p  -  p  =  .10 
1    2 

r  -  2/hour 

1 

r  =  1/hour 
2 

Table  8  lists  the  results  for  several  possible  values  of  V  *. 

r 

While  the  scope  of  these  investigations  of  value  sensitivity  is 
obviously  limited,  it  appears  that  over-estimation  of  the  value  of 
targets  of  higher  worth  might  impose  less  risk  than  under -estimation. 
Further  examination  in  the  area  of  sensitivity  where  several  weapon 
types  are  employed  simultaneously  would  be  necessary  before  conclusive 
results  could  be  stated. 
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If  V  *  is  misestimated  by  a   factor  of   z, 

v  * 

r 

the 
z=10 

y 

expected   returns  decrease  by  y. 
z=5              z=l/2              z=l/10 

y             y                y 

1 

287. 

227. 

57. 

87. 

2 

8% 

67. 

107. 

377. 

3 

47. 

37. 

57. 

537. 

5 

27. 

17. 

37. 

677. 

7 

17. 

17. 

17. 

667. 

10 

17.             <.57. 

• 

17. 

407. 

TABLE   8 
PERCENTAGE  DECREASE  IN  EXPECTED  RETURNS 
DUE  TO  MIS  ESTIMATION  OF  V  * 
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4.   CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  for  determining  an  optimal  initial  loading  for  missions 
against  targets  of  opportunity  has  been  developed.  As  a  by-product 
of  this  load  determination,  an  optimal  strategy  for  the  sequential 
expenditure  of  weapons  has  been  found.   The  gains  in  the  expected 
value  of  targets  killed,  when  the  initial  loading  is  optimal  and  when 
an  optimal  strategy  is  employed,  have  been  shown  to  be  significant 
when  compared  with  initial  loadings  and  firing  doctrines  derived  by 
the  application  of  simple  intuition  and  judgement  to  kill  probabilities, 
target  values,  and  encounter  rates. 

The  determination  of  an  optimal  loading  and  expenditure  policy, 
however,  are  at  present  of  such  complexity  as  to  require  the  use  of  a 
high  speed  computer  in  all  but  the  most  trivial  cases.  For  those 
operational  situations  in  which  a  computer  will  not  be  available,  a 
hand-calculable  approximation  method  for  the  determination  of  a  best 
initial  loading  has  also  been  developed,  which,  though  sub-optimal, 
has  been  found  to  generate  loads  which  yield  high  expected  returns. 

Either  of  the  above  methods  should  yield  significantly  better 
results  than  present  policies  which  prescribe  that  aircraft  be  loaded 
for  targets  which  are  assumed  known  or  that  several  aircraft  each  be 
given  different  loads  in  order  to  cover  an  unknown  spectrum  of  targets. 

The  models  presented  in  this  paper  assume  that  weapons  can  be 
loaded  in  any  combination  which  satisfies  rack  and  weight  constraints, 
and  that  the  pilot  has  complete  freedom  of  choice  as  to  the  sequence 
in  which  weapons  can  be  expended.  For  most  attack  aircraft  these 
assumptions  are  not  always  true.  Restrictions  on  unsymmetrical  loadings 
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are  one  limiting  factor.   Others  relate  to  the  mechanics  of  multiple 
bomb  racks.  The  models  should  be  expanded  to  account  for  these  aspects 
of  the  problem. 

A  vital  assumption  in  the  developments  of  both  models  was  that 
inter-encounter  times  would  be  distributed  exponentially.   It  is 
recommended  that  an  analysis  of  operational  data  be  made  to  ascertain 
the  validity  of  this  supposition. 

It  seems  reasonable  to  assume  that  a  list  of  all  possible  types 
of  targets  of  opportunity,  for  a  spectrum  of  scenarios,  would  be 
small  enough  -  certainly  not  larger  than  50  -  that  a  thorough  analysis, 
as  suggested  in  Section  3,  could  be  made  of  their  relative  values.   Such 
an  analysis  would  indicate  general  strategic  values  which  could  then  be 
adjusted  for  specific  tactical  scenarios.  Furthermore,  additional 
investigations  similar  to  that  outlined  in  Section  3  should  be  under- 
taken to  study  the  sensitivity  of  the  optimal  load  and  strategy  cal- 
culations to  these  values. 

The  hand  approximation  method  is  rather  crude  in  form.   It  is 
recommended  that  this  model  be  examined  further  in  an  attempt  to  refine 
the  calculations  and  to  extend  the  model  to  cover  a  wider  variety  of 
cases  in  a  more  complete  fashion. 

The  computer  programs  of  Appendix  II  were  written  solely  for  the 
computations  used  in  this  thesis  and  should  be  improved  upon  before 
extensive  investigations  are  made  of  the  dynamic  programming  treatment 
of  the  problem.  The  four  target-seven  weapon  example  required  two 
hours  and  fifteen  minutes  on  a  CDC-1604  computer  for  the  determination 
of  an  optimal  load.  This  is  testimony  to  the  inefficiency  of  the 
"BESTLD"  FORTRAN  program. 
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The  problem  of  determining  an  optimal  initial  loading  is  similar 
to  the  classic  "flyaway-kit"  problem  of  dynamic  programming.   In  the 
strike  aircraft  loading  problem  the  additional  complication  is  present 
that  the  costs  of  shortages  are  not  constant  but  depend  on  the  initial 
loading.   It  would  seem  though  that  research  into  solution  of  the 
strike  aircraft  problem  by  this  method  should  be  fruitful.   Very  good 
approximations  to  an  optimal  initial  loading  should  be  obtainable  very 
rapidly  in  this  way. 

Although  the  "POLICY"  FORTRAN  program  computes  the  optimal 
strategy  with  relative  ease,  usually  less  than  four  minutes  is  required 
on  a  CDC-1604  computer,  it  is  recommended  that  methods  of  approximation 
be  developed  for  which  a  strategy  can  be  calculated  without  the  aid  of 
a  high  speed  computer. 


1 

Bellman,  R.  E.  and  Dreyfus,  R.  E.,  Applied  Dynamic  Programming, 
(Princeton,  New  Jersey:   Princeton  University  Press,  1962),  pp.  42-43. 

47 


ACKNOWLEDGEMENTS 

The  authors  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  Stephen  M.  Pollock  of  the 
U.  S.  Naval  Postgraduate  School  and  to  Professor  David  A.  Schrady 
whose  suggestions  and  criticisms  helped  determine  the  form  of  the 
final  write-up. 


48 


BIBLIOGRAPHY 

1.  Bellman,  R.  E.   Dynamic  Programming.  Princeton  University  Press, 

1957. 

2.  Bellman,   R.    E.    and  Dreyfus,    S.    E.     Applied  Dynamic   Programming. 

Princeton  University  Press,    1962. 

3.  Parzen,    E.      Modern  Probability  Theory  and  Its  Applications. 

John  Wiley  and  Sons,    1960. 

4.  Stevens,    S.    S.      A  Metric    for   Social  Concensus,    Science,    vol.    151, 

4  February   1966,    pp.    530-541. 

5.  Bram,    J.      Allocation  of  Weapons   to  Targets  with   Exponential 

Arrival  Times   in  a   Limited  Time   Interval.      Operations  Evaluation 
Group,    Center   for  Naval  Analysis.      Interim  Research  Memorandum 
32,    1963. 

6.  U.    S.    Naval  Ordnance  Test  Station.      Preferred  Aircraft  Loadings 

Based  on  the  Predicted  Combat  Effectiveness  of  Conventional 
Air-to-Surface  Weapons.  NAVWEPS  Report  8664,  NOTS  TP  3710, 
April   1965.    SECRET 


49 


APPENDIX  I 

ENCOUNTER  RATES 

In  a  combat  situation,  many  targets  (especially  "targets  of 
opportunity")  may  be  thought  of  as  being  distributed  randomly  through- 
out a  given  area,  in  that  the  pilot  flying  the  mission  does  not  know, 
beforehand,  where  or  when  a  target  will  be  encountered.  This  may  be 
thought  of  as  the  classical  situation  in  which  it  may  be  assumed  that, 
in  any  small  time  interval  dt,  the  probability  that  a  target  is  en- 
countered is  proportional  to  dt  and  is  independent  of  any  previous 
encounters. 

Over  a  period  of  time,  photo  or  visual  reconnaissance  can  give 
an  average  density  of  targets  in  a  given  area.  From  the  attack  air- 
craft speed  and  the  average  search  sweep  width  (width  of  a  path  in 
which  the  pilot  could  be  expected  to  see  a  target)  the  area  swept  per 
unit  time  can  be  obtained.  By  multiplying  the  average  type  i  target 
density  by  the  average  sweep  area  per  unit  time,  an  estimate  of  the 
average  encounter  rate  r.  is  obtained  for  target  type  i. 

The  average  encounter  rate  developed  above  can  be  used  to  gen- 
erate the  distribution  of  inter-encounter  times. 

Assumptions: 

Probability  of  one  encounter  in  dt  =  r.dt  +  h(dt)  where  h(dt) 
approaches  zero  as  dt  approaches  zero. 

Probability  of  two  or  more  encounters  in  dt  »  h(dt) 

Probability  of  no  encounters  in  dt  ■  1-r.dt  +  h(dt)  p(t)dt  ■ 

prob.  (t  ^  first  encounter  time  ^  t  +  dt)  =  dP(t) 

P(t)  =   P(T)dt  =  prob. (encounter  time  ^  t) 
o 
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Q(t)  ■=  l-P(t) 

Thus  tf  t   is  the  encounter  time,  and  t  is  any  arbitrary  time  in 

the  mission, 

Q(t+dt)  -  prob.((t  >  t)(no  encounter  in  dt)) 

e 

-  Q(t)  (1-r.dt  +  h(dt)) 

Subtracting  Q(t)  from  both  sides,  dividing  through  by  dt  and  taking 

the  limit  as  dt  approaches  zero  gives  dQ(t)/dt  =  -Q(T)r 

Rearranging  gives  dQ(t)/Q(t)  ■  -r  dt 

i 

Integrate  both  sides  to  get  In  Q(t)  =  C-r  t 
Thus  Q(t)  =•  Ae"*!* 

But  Q(0)  ■  1  since  there  can  certainly  be  no  encounters  in  zero  search 

-r.t  -r.t 

time.  Therefore  A  ■  1  and  Q(t)  =  e   x  and  consequently  P(t)  =  l-e  *■ 

~r  t 
and  p(t)  ■  re   i 

Thus  it  is  seen  that  the  inter-encounter  times  are  exponentially 

distributed.   It  also  follows  naturally  that  the  probability  distribution 

on  the  number  of  targets  type  i  encountered  during  a  mission  of  length 

T  is  Poisson,   i.e. , 

n  -r  .T 
prob.  (N(T)  =  n)  =  (r.T)   e  x   /n! 

and  the  expected  no.  N(T)  =  r  T   (This  is  obtained  from  the  observations). 

i 


Parzen,  E.,  Modern  Probability  Theory  and  Its  Applications, 
(New  York:  John  Wiley  and  Sons,  I960)  p.  253. 
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APPENDIX  III 
TABULATED  VALUES  OP  f 


ijs 


BJS   RT  *  .05 


10 


.15 


.20 


1 

.C4877 

.09516 

2 

.04998 

.09984 

3 

.05000 
.05000 

.10000 

4 

.10000 

5 

.05000 

.10000 

6 

.05000 

.10000 

7 

.05000 

.10000 

8 

.05000  . 

.  10000 

9 

.05000 

.10000 

10 

.0  5000 

.10000 

11 

.05000 

.10000 

12 

.05000 

.10000 

13 

.05000 

. 10000 

14 

.0  5000 

.10000 

15 

.05000 

.10000 

16 

.05000 

.10000 

17 

.05000 

.10000 

18 

.05000 

.10000 

19 

.05000 

.10000 

20 

.05000 

.10000 

13929 

.18127 

.22120 

14948 

.19879 

.24770 

14998 

.  19994 

.24986 

15000 

.20000 

.24999 

1  5000 

.20GOO 

.25000 

15C00 

.20000 

.25000 

15000 

.20000 

.25000 

15000 

.20000 

.2  5000 

15000 

.20000 

i. 25000 

15000 

.20000 

1  .25000 

15000 

.20000 

.25000 

15000 

.20000 

.25000 

15000 

.20000 

.25000 

15000 

.20000 

..25000 

15000 

.20000 

.25000 

15000 

.20000 

.25000 

15000 

.20000 

.25000 

15000 

.20000 

.25000, 

15000 

.20000 

.25000 

15000 

.20000 

.25000 

BJS   RT  =  .30 


.35 


.40 


1 

.25918 

.29531 

.52968 

2 

.29612 

.34398 

.39123 

3 

.29972 

.34949 

.39916 

4 

.29998 

.34997 

.39993 

5 

.50000 

.35000 

.40000 

6 

.30000 

.35000 

.40000 

7 

.30000 

.35000 

.40000 

8 

.30000 

.3500  0 

.40000 

9 

.30000 

.35000 

.40000 

10 

.30000 

.35000 

.40000 

11 

.30000 

.35000 

.40000 

12 

.30000 

.35000 

.40000 

13 

.30000 

.35000 

.40000 

14 

.30000 

.35000 

.40000 

15 

.30000 

.35000 

.40000 

16 

.30000 

.35000 

.40000 

17 

.30000 

.35000 

.40000 

18 

.30000 

.35000 

.40000 

19 

.30000 

.35000 

.40000 

20 

.30000 

.35000 

.40000 

.45 

.36237 

.43781 

.44869 

.44989 

.44999 

.45000 

.45000 

.4  5000 

.4  5000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45000 

.45C00 

.45000 


.50 

.39347 
.48567 
.49806 
.49981 
.49998 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 
.50000 


TABLE  11 
VALUES  OP  f 


ijs 
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BJS   RT  =  .55 


1 

2 
3 
M 
5 
6 
7 

B 
9 
10 
1  1 
12 
13 
14 
15 
16 
17 
18 
19 
20 


.42305 
.52  078 
.54724 
.54971 
.54997 
.55000 
•55000 
.5  5000 
.5  500  0 
.55000 
.55000 
.55000 
.55^00 
•55000 
.5500  0 
.55000 
.55u00 
.55000 
•55000 
.55000 


.60 


.65 


.451 19 

.47795 

.57  30  9 

.61658 

.59621 

.64492 

.59956 

.64937 

.59996 

.64993 

.600  00 

.64999 

.6000  0 

.65000 

.60000 

.65000 

.60000 

.6  500  0 

.60000 

.65000 

.60000 

.6  5000 

.60000 

.65000 

.6000u 

.6  5000 

.60000 

.65000 

.60000 

.65000 

.60000 

.65000 

.60000 

.65000 

.60000 

.65000 

.60000 

.6  5000 

•60n00 

.65000 

.70 

50341 
65922 
69336 
6991  1 
6999o 
69999 
70(  00 
70000 
70000 
7000  0 
70000 
70u00 
70000 
70000 
70000 
70000 
70000 
70000 
70000 
70C00 


.75 

.52763 
.70099 
.74150 
.74879 
.74985 
.74998 
.75000 
.75000 
(.75000 
'.75000 
.75000 
.75000 
.75000 
.75000 
.75000 
.75C0O 
.75000 
.75000 
.75000 
.75000 


BJS   RT  =  .80 


.65 


.90 


.95 


1 

.55067 

.57159 

.59343 

.6  1326 

2 

.74188 

.78187 

.62095 

.8591 1 

3 

.78930 

.83675 

.68580 

.93045 

4 

.79638 

.84788 

.69726 

.94653 

5 

.79979 

.84971 

.69961 

.94948 

6 

.79998 

.84997 

.89995 

.94995 

7 

.8000  0 

.8  5000 

.69999 

.94999 

8 

.80000 

.8  5000 

•90000 

.95000 

9 

.80000 

.85000 

.90000 

.95000 

10 

.80000 

.8  5000 

'  .90000 

.95000 

11 

.80000 

.85000 

.90000 

.95000 

12 

.80000 

.8  50 00 

.90000 

.95000 

13 

.80000 

.85000 

.90000 

.95000 

14 

.80000 

.85000 

.90000 

.95'  00 

15 

.80000 

•85C00 

.90000 

.95000 

16 

.80000 

.85000 

.90000 

.95000 

17 

.80000 

.85000 

.90000 

.95000 

18 

.80000 

.85000 

.90000 

.9  5000 

19 

.80000 

.65000 

.90000 

.95000 

20 

.80000 

.85000 

.90000 

.95000 
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BJS 


Rl 


2 


1 

2 

3 

4 

5 

6 

7 

6 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


632 
896 
976 
995 
999 
999 
999 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 


12 
36 
66 
65 
31 
91 
99 

ou 

00 
00 
00 
00 
00 
00 
00 
00 
00 

ou 

00 
00 


86460 
45866 
78198 
92486 
97751 
994o6 
99861 
99971 
99994 
99999 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 


95021 
75106 
32787 
68064 
86538 
94930 
98281 
99471 
99851 
99962 
99991 
99998 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 


98168 
89011 
65200 
21853 
58970 
80457 
91524 
96637 
98774 
99587 
99871 
99962 
99990 
99997 
99999 
00000 
00000 
00000 
00000 
00000 


99326 
95283 
82818 
56316 
12266 
5067U 
74452 
87769 
94598 
97761 
99151 
99696 
99898 
99968 
99990 
99997 
99999 
00000 
00000 
00000 


8JS 

1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
18 
17 
18 
19 
20 


RT  = 


1 
2 
3 
4 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
o 


,99752 
,9801? 
91820 
76700 
48194 
03626 
42996 
68596 
83874 
92267 
96529 
98538 
99421 
99763 
99923 
99974 
99992 
99998 
99999 
00000 


99909 
99179 
96216 
86039 
70740 
40669 
95698 
35827 
62918 
79866 

89720 
6.95o55 
6.97755 
6.99o36 
6.996o6 
6.99846 

99944 

99980 
99993 
99998 


8 

99966 
99665 
98269 
94051 
84068 
64964 
33627 
88331 
29076 
57414 
7582? 
87017 
93397 
96815 
96541 
99364 
99736 
99895 

99960 
99966 


99988 
99864 
99241 
97118 
91622 
5.80053 
6.59375 
26965 
81420 
22679 
5206U 
71779 
84202 
91587 
95734 
97937 
99048 
99580 
99823 
99928 


10 

99995 
99946 
99669 
98635 
95710 
89001 
75987 
53965 
8.20663 

8.74890 
9.16586 
46908 
67753 
81306 
89652 
94526 
97230 
98658 
99377 


9.99722 
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BJS 


RT  = 


11 


12 


13 


14 


15 


l 
2 
5 
4 

5 
6 

/ 

B 

V 

10 

1 1 

12 
13 
14 
15 
16 
W 
18 
19 
20 


9999b 
9  9  9/6 


2i9985/ 
3. ,99366 
4C97855 


5 
6 

/ 

6 
9 
9 
10 
10 
10 
10 
10 
10 
10 
10 
10 


941Q3 
86242 
71923 
48724 
1  4  6  7  v5 
68664 
10757 
41866 
63/3* 
78334 
87615 
93207 
9  6  4  2  6 
98195 
9912^ 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

10 

11 

11 

11 

11 

11 

11 

11 

11 


99999 
99991 
9993V 
9971U 
96950 
96916 
92334 
83383 
67860 
43641 
0891b 
62759 
05162 
37006 
59806 
75364 
65493 
91/90 
95532 
97660 


1 
1 
2 
3 
A 
5 

6 

7 

8, 

9, 

10. 

11, 

11, 

11  . 

12. 

12. 

12. 

12. 

12. 

12. 


•00000 
,9999/ 
,999/4 
99869 
99495 
98422 
95633 

90430 
8045b 
63873 
38/05 
0^369 
5  70/8 
9  9  7  7  4 
32261 
55900 
72331 
83304 
9028/ 
94534 


1 
1 
2 
3 

4 

5 
6 

7 

8 
9 
10 
11 
11 
12 
12 
13 
13 
13 
13, 


00000 
99999 
99969 
99942 
99761 
99208 
97765 
94623 
88418 
77476 
59910 
33906 
98060 
5161b 
94572 
27636 
52044 
6932^ 
81059 


13.68710 


1 
1 
2 
3 

5 

7 

8 

9 
10 
11 
12 
12 
13 
13 
14 
14 
14 
14 


00000 
99999 
99996 
99974 
99869 
99610 
98846 
97046 
93301 
86316 
74470 
55994 
29233 
92912 
46346 
89537 
23125 
48239 
66292 
78770 


BJS 


RT  =  16 


17 


18 


19 


20 


1 
2 
3 
4 
5 
6 
7 

b 

10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 


1 
2 
2 
3 

o 
7 

8 

9 
10 
11 
12 
13 
13 
14 
14 
15 
15 


00000 

oooou 

9999b 
99969 
99949 
99810 
99410 
9841u 
96211 
91681 
84142 
71442 
52130 
24679 
87926 
41252 
84656 
18721 
44486 


1 

d 

I 

3 

4 

5 

o 

7 

6 

9 

10 

11 

12 

13 

14 

14 

15 

15 

16 


00000 
00000 
99999 
99995 
9  9  9  7  7 
99909 
99703 
99160 
97900 
95286 
90373 
81906 
68406 
48319 
20235 
83090 
36316 
79914 
14416 


15.63262   16.40786 


1 
2 
5 
3 
4 

6 

7 

6 

9 

10 

11 

12 
13 
14 
15 
15 
16 
16 
17 


•0J0U0 
.  00000 
•00000 
.  99996 
, 99990 
.9995/ 
.99853 
.9956*5 
.98858 

.97320 
.94283 
.88794 
.79627 
.65368 
.44560 
.15895 
.78390 
.31525 
.75300 
.10209 


1 
2 
3 
3 
4 
5 
6 
7 

8 
9 

10 
11 
12 
13 
14 
15 
16 


00000 
00000 
00000 
99999 
99995 
99980 
99928 
9  9  7// 
99369 
98504 
9  6  6/2 

93204 
87146 
7  7  3  o  6 
62333 
4o854 
11651 
16.73814 
17.26866 
17.70805 


1 
2 

3 

4 
4 

3 

b 

7 

8, 

9, 

10. 

11, 

12. 

13. 

14. 

15. 

16. 

17. 

17. 


•00000 

•00000 

-00000 

.00000 

,99998 

99991 

99965 

99887 

99679 

99179 

98098 

95959 

92058 

85445 

74959 

59306 

37200 

07497 

69355 


18.22329 


67 


INITIAL  DISTRIBUTION  LIST 


No.   Copies 


1.  Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  20 

2.  Library 

U.  S.  Naval  Postgraduate  School 

Monterey,  California  2 

3.  Chief  of  Naval  Operations 

Office  of  Naval  Warfare  Analysis  (OP-91) 

Department  of  the  Navy 

Washington,  D.  C.  1 

4.  Captain  J.  Cummings 
Office  of  Program  Appraisal 
Department  of  the  Navy 

Washington,  D*  C.  1 

5.  J.  Bram 

Operations  Evaluation  Group 

Office  of  the  Chief  of  Naval  Operations 

Washington,  D.  C.  1 

6.  Carl  L.  Schaniel 

Naval  Ordnance  Test  Station 

China  Lake,  California  1 

7.  Prof.  Stephen  M.  Pollock 
Department  of  Operations  Analysis 
U.  S.  Naval  Postgraduate  School 

Monterey,  California  2 

8.  Prof.  David  Schrady 
Department  of  Operations  Analysis 
U.  S.  Naval  Postgraduate  School 

Monterey,  California  1 

9.  Prof.  Roger  Schroeder 
Department  of  Operations  Analysis 
U.  S.  Naval  Postgraduate  School 

Monterey,  California  1 

10.  LT  Thomas  M.  Burns 
62  Work  Avenue 

Del  Rey  Oaks,  California  1 

11.  LT  John  N.  Donis 

230  Grove  Acres,  Apt.  332 

Pacific  Grove,  California  1 

68 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  claealllcatlon  ol  title,   body  ol  abatract  and  Indexing  annotation  mutt  be  entered  when  the  overall  report  la  clmaalllad) 


1     ORIGINATING  ACTIVITY  (Corporate  author) 

U,    S.    NAVAL   POSTGRADUATE   SCHOOL 
MONTEREY,    CALIFORNIA 


2a      REPORT    SECURITY     CLASSIFICATION 

UNCLASSIFIED 


2b     CROUP 


3     ntPOHT    TITLE 


OPTIMUM  AIRCRAFT  WEAPONS  SELECTION  AND  SEQUENTIAL  FIRING  STRATEGY  FOR 
TARGETS  OF  OPPORTUNITY 


4     DESCRIPTIVE   NOTES  (Type  ol  report  and  Incluelve  datee) 

Master's  Thesis 


S     AUTHORS)  (Laet  name,    llrat  name,    Initial) 

Donis,  John  N.  and  Burns,  Thomas  M. 


6     REPORT   DATE 


May    1S66 


la       TOTAL.    NO      OF     PACES 

67 


7b.     NO      OF    REFS 


8a      CONTRACT    OR    GRANT    NO. 


b      PROJEC  T    NO 


9a     ORIGINATOR'S    REPORT   NUMBERCS; 


9b     OTHER  REPORT    noCSJ    (Any  other  numbera  that  may  be  aaalgned 
thla   report) 


10     AVAILABILITY/LIMITATION   NOTICES 


11     SUPPLEMENTARY    NOTES 


12     SPONSORING   MILITARY  ACTIVITY 


U.  S.  Naval  Postgraduate  School,  Monterey, 
California 


13     ABSTRACT 


A   pilot   engaged    mi  a   strike  mission  against   targets   of   opportunity   is   cu   - 
carnally  confronted  with   the  decision  as   to  he  j      n     whicl    o  s   available 

weapon;    should   bo   cxp  iiidcd  u   u   v.  riety   ui    targets  which   are   encountered. 

Associated   ,/il.u    this  decision  problem   is    the  problem  of   determining   an  optimal 
loud   lor   such  a   mission  subject   to   the   constraints   of  available   payload  and 
number  of  weapons   stations  on  the  aircraft.      Certain  assumptions  are  made  con- 
cerning  the  distribution  of   targets  within   the   target   are:   which    lead    to   a 
dynamic   programming   formulation  of   the   decision   problem.      This   yields   a    system 
of   ordinary  differential   equations  which  are   solvable   recursively.      In  addition 
to    the  dynam-c   programming  model,    a   sub-optimum  determination  of   a     'best    load" 
is  discussed.      Although   nut   a^    c<-..i    lete   or      recise  as   the  dynamic   programming 
method,    this    formulation  .  1.  iptcl  lo    to   squadron   level  decision 

making. 
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